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Abstract: We derive the evolution equation for hadronic scattering amplitude at high energy. Our 
derivation includes the nonlinear effects of finite partonic density in the hadronic wave function as 
well as the effect of multiple scatterings for scattering on dense hadronic target. It thus includes 
Pomeron loops. It is based on the evolution of the hadronic wave function derived in (T|. The kernel 
of the evolution equation defines the second quantized Hamiltonian of the QCD Reggeon Field 
Theory, Hrft beyond the limits considered so far. The two previously known limits of the evolution: 
dilute target (JIMWLK limit) and dilute projectile (KLWMIJ limit) are recovered directly from 
our final result. The Hamiltonian Hrft is applicable for the evolution of scattering amplitude for 
arbitrarily dense hadronic projectiles/targets - from "dipole-dipole" to " nucleus- nucleus" scattering 
processes. 
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1. Introduction 

This paper continues our investigation of high energy evolution of hadronic cross sections and other 
observables. The subject has a long history, starting with Gribov's ideas of reggeon field theory ||] 
and early work exploring Pomeron interactions Q, [Q]. In the framework of QCD the perturbative 
evolution equation - the BFKL equation was derived in |5). It describes the evolution of forward and 
non-forward scattering amplitudes as the energy of the collision becomes large. The BFKL equation 
was clearly a crucial milestone in the study of high energy scattering. It has given an impetus to a 
lot of theoretical and experimental work. However it has been realized very quickly that at very high 
energies the BFKL equation violates unitarity and leads to cross sections which rise as a power of 
energy. This violation is very severe, since it does not only violate the Froissart bound for the total 
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cross section, but also violates unitarity of the scattering amplitude at fixed impact parameter. In 
order to restore unitarity one has to include i-channel exchanges with more than just two reggeized 
gluons. This idea has been developed by Bartels |J. Within this framework the generalization of 
the BFKL equation is the BKP equation, which governs (perturbative) high energy evolution of an 
amplitude due to an exchange of an arbitrary but fixed number of reggeized gluons (7j. 

It was further realized that to achieve the s-channel unitarity one has to introduce transitions 
between states with different number of reggeized gluons ||. The ideas put forward in [|j have 
been under intense investigation during the last decade or so |, H, 0, U, H, Q. At present, this 
approach provides elements of an effective theory in terms of t-channel gluon states and transition 
vertices. These elements have been put together in [14] into an effective theory of BFKL Pomerons 
interacting via triple pomeron vertex of ||. This model is meant to describe nucleus-nucleus collisions 
at high energies at LLA and large N c . We also note that Lipatov and collaborators have 
derived an effective action with both real and reggeized gluons as effective degrees of freedom. This 
action respects the unitarity of full QCD, but its complexity has so far precluded any progress in 
understanding its physical consequences. 

It is hoped that the end result of this direction of research will be a quantum theory of interacting 
Reggeons derived entirely from QCD. 

A parallel line of research originated by Gribov, Levin and Ryskin |]l6| has been vigorously 
pursued in the last 15 years or so. It is based on the idea of gluon saturation. At high energy the 
evolution of physical observables should slow down because of nonlinear effects due to large density of 
gluons in the hadronic wave function. The GLR equation - the nonlinear evolution equation for gluon 
density in the double logarithmic approximation was derived in |l(| and put on firmer theoretical 
grounds in [17]. The gluon saturation ideas have been further developed in a series of papers by 
Mueller [18], who also introduced the notion of QCD dipoles as a convenient basis for the discussion 
of high energy processes, and has related it to BFKL Pomeron and the triple Pomeron vertex [19]. 

It was noted in |2(J that the problem of saturation can be also related to nonlinearities in the 
classical Yang-Mills equations. Following this observation, the formal path integral approach to 
the problem of evolution has been developed in pi]]. This together with independent approaches 
of |2^] and |2^| , and a later work [^4|] resulted in the derivation of a functional evolution equation 
for the correlators of the color charge density in the hadronic wave function, the so called JIMWLK 
equation (also called sometimes the Balitsky hierarchy) . The equations are more conveniently written 
as evolution of the correlators of Wilson lines, introduced in [25] as effective degrees of freedom at high 
energy. The JIMWLK equation takes into account the effects of nonlinearities in the wave function 
of the projectile, but when applied to the scattering amplitude, it does not account for the bulk 
of multiple scattering corrections. It is applicable in a situation when a dense object ("nucleus") 
scatters off a dilute perturbative object ("dipole"). The equation therefore does not include the 
Pomeron loop effects [26, 27, |2^, pOfl . The reverse situation - that of a perturbative "dipole" 
scattering off a dense "nucleus" was considered in [31| and the evolution equation for this limit has 



been derived. We will refer to this as the KLWMIJ equation [31]. Although it describes the same 



physical process as JIMWLK, the KLWMIJ evolution acts on the wave function of a dilute object, 
and thus does not include the effects of nonlinearities. It does, on the other hand, take into account 
all important multiple scattering effects. The efforts to consistently include Pomeron loops into the 



- 2 - 



evolution have continued since, yielding many interesting developments [32, |33|, [?4|, 35, 36, 37 



39, 40, 41, 42, [43, 0JJ]. Also a statistical analogy of the Pomeron loop effects has been suggested in 



|45, 46], although its validity has not been convincingly demonstrated. 

Other recent interesting developments in this area include the calculation of the next to leading 



order corrections to the JIMWLK kernel p7j, 48] as well as the generalization of the approach to 



calculate less inclusive observables than the scattering amplitude ]^9|, |5(], |5T], ^] . 

The last years also have seen progress in relating the JIMWLK approach to the reggeon field 
theory. This relation has been discussed in some detail in [j53|| . The functional JIMWLK equation 
has direct interpretation as the Schroedinger equation of the Hamiltonian reggeon field theory. The 
JIMWLK kernel in this view is simply the Hamiltonian of the RFT (Reggeon field theory). To 
be precise, the theory which we dub RFT is the effective theory of QCD at high energies which 
governs the high energy behavior of scattering amplitudes of hadrons. The degrees of freedom of 
this effective theory are the Wilson lines (single gluon scattering amplitudes). These are not quite 
the reggeized gluons of f|],[[?J but are closely related to them. The mapping between the Wilson 
lines and reggeized gluons of the standard formulation has been discussed in detail in [53]. Since 



JIMWLK equation does not contain the Pomeron loop effects, the RFT Hamiltonian obtained in this 
way is not complete. To obtain the complete RFT Hamiltonian one has to take into account both the 
effects of the nonlinear corrections in the projectile wave function as well as contributions of multiple 
scattering corrections to the scattering amplitude, in the same framework. Although the complete 
RFT Hamiltonian has not been derived yet, some of its general properties have been discussed in 



1 54]. In particular, since Hrft generates the energy evolution of hadronic scattering amplitudes, 



unitarity requires the spectrum of its eigenvalues to be positive. 

An important step towards the derivation of the complete Hrft was made in Q. This paper 
derived the evolution of the hadronic wave function under boost. The one step that is still missing 
in H is the derivation of Hrft itself. To derive Hrft one has to consider calculation of physical 
observables, such as the scattering amplitude in the evolved wave function, rather than simply the 
evolution of the wave function itself. 

The main purpose of the present paper is precisely to make this additional step and complete 
the derivation of Hrft- Our final result for the RFT Hamiltonian is given in eq. (pL2Sl) . This result 
is valid to leading order in the coupling constant a s at any parametric value of the color charge 
density of the colliding hadrons. It also accounts for all multiple scattering effects in the eikonal 
approximation. The term eikonal is used here only in the sense that individual partons in the 
hadronic wave function scatter eikonally. The hadron as a whole in this approximation does undergo 
both elastic and inelastic scattering fl55[| . This approximation is the same as used in all the recent 
papers on the subject quoted above. 

The Hamiltonian Hrft effectively resums all perturbative diagrams in which every factor a s 
is enhanced by a single logarithm of energy (leading log approximation). In all regimes explored 
previously, in addition to the leading log approximation some assumption about partonic densities was 
made. The density is either considered to be parametrically small (dilute regime - KLWMIJ evolution) 
or parametrically large (dense regime - JIMWLK evolution) . Although we will be frequently referring 
to the dense regime in this paper, the Hamiltonian we derive is valid for arbitrary density and thus 
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has only one resummation parameter a s Y ~ 1. Corrections to our result are suppressed by a power 
of a s at any partonic density. Although we have not studied these corrections in full generality, one 
can easily see that in both the JIMWLK and the KLWMIJ limits the expansion is in integer powers 
of a s and any correction to our result contributes at most at next-to-leading order in this expansion. 



The result presented in eq.( 4.28 ) generalizes all previously available limiting expressions for Hrft 
HI, ||, ||, ||, H, ||, ||, || to the most general situation. It reproduces the JIMWLK and the 
KLWMIJ Hamiltonians in the appropriate limits, and is applicable also when both colliding objects 
are either small ("dipole-dipole scattering") or large (" nucleus- nucleus scattering"). 



Similarly to previously known limits, the complete Hrft eq.( 4.28 ) can be expressed in terms of 
Wilson lines. Those are SU(N) unitary matrices which depend on the two dimensional transverse 
coordinates. These are the basics quantum fields - the degrees of freedom of the effective Reggeon 
Field Theory. We note that one can introduce two sets of Wilson lines (associated to the propagation 
of projectile or target partons), which play the role of (almost) canonically conjugate variables. Either 
set of the Wilson lines can be chosen as a complete set of variables on the Hilbert space of RFT. 

The scattering amplitude is not the only physical observable whose evolution is governed by 
Hrft- Ln the companion paper [56| we continue this line of investigation by deriving expressions for 
more exclusive observables, which also resum effects of Pomeron loops. In particular we discuss in 
1 56] the single gluon inclusive spectrum and also the double and multi gluon exclusive spectra when 
none of the observed gluons are separated by large rapidity, so that the effects of energy evolution 
between the final state gluon rapidities are unimportant. 

This paper is structured as follows. In Section 2 we recap the main results of Q and outline the 
steps for calculation of Hrft- Section 3 contains the calculations of main ingredients necessary for 
the final push. In Section 4 we complete the calculation of the Hrft and derive our main result. 
We discuss our results in Section 5. Appendices provide some details of calculation as well as discuss 
some aspects of the perturbative expansion of our result. 



2. The Summary and the Road Map. 

In we have analyzed the evolution of the hadronic wave function at high energy. Here we summarize 
the results of Q. 

The standard expansion of the field operator in terms of creation and annihilation operators of 
free gluons is 

f°° dh^ rl 2 k 1 r _l . -i 

At{x-,x) = J ^^-^={aUk + ,k)e-* k+ *--^ + a? (k+ , k) e* k+ *~ +* k *} (2.1) 
where 

[a?(k + ,k),af(p + ,p)] = (2Trf5(k + -p + )5 2 (k-p). (2.2) 

The operator A°^(x~,x) is the vector potential operator in the light-cone gauge A + = with the 
zero longitudinal momentum mode subtracted. 1 Throughout the text we will be using upper Latin 

it was denoted as A. 
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indices a,b . . . = 1, N 2 — 1 to denote color while lower indices = 1, 2 to denote rotational 

components. For further reference we note that the vector potential can also be decomposed in the 
creation and annihilation operators in the frequency basis: 

I" 00 rlk~ rl 2 k 1 ( ■ k 2 ; ■ k 2 , } 

An/ — \ I u <" I a/i — i \ —i— — x —ikx , at/r — n i— — x +ikx I / n n \ 

A i( x > x ) = I 4^2 ^||= ' fc ) e 2fc " + a i ( fc ,fc)e»- | (2.3) 



o 



where 



so that 



<{k-,k) = yj!£-a$(k + ,k) (2.4) 

[at(k-,k),af(p-,p)} = (2nf5(k--p~)5 2 (k-p). (2.5) 

We consider a hadronic projectile moving to the right with large energy. Suppose at initial 
rapidity we know the hadronic wave function The bulk of gluons in this wave function are at 
rapidities greater than some "cut off' rapidity. The properties of this valence component of the wave 
function are characterized by correlators of the color charge density operator j a (x). The gluons with 
rapidity smaller than the cutoff are small in number and therefore do not contribute to observables 
in the leading order in the coupling constant. 

This "soft" component of the wave function can be calculated perturbatively. It was shown in 
HI that the wave function, including its soft component has the form 

|*) = n[a,<J,j]\v) (2.6) 

where a and are soft gluon creation and annihilation operators with rapidities below the cutoff. 
The valence state \v) has no soft gluons and is therefore annihilated by the soft gluon annihilation 
operators 

a\v) =0. 

The evolution operator O is a unitary operator of the Bogoliubov type 

fl = CB. (2.7) 
Here C is a coherent operator that creates the "classical" Weizsaker- Williams field 

C = exp|2i J d 2 xb a i {x)A a i {x- = 0,x)| . (2.8) 

and B is a Bogoliubov type operator responsible for the leading quantum corrections (see below). 

The color charge operators j(x) are operators on the valence Hilbert space and form a local 
SU{N) algebra: 

[f(x),j b (y)] = i 9 r bc f(x)5 2 (x-y). (2.9) 

The Weizsaker- Williams field bf(x) created by the valence modes depends only on transverse coordi- 
nates x and is a two dimensional pure gauge field. It can therefore be written in terms of a unitary 
matrix U as 
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where f abc are the structure constants of the SU(N) and U ab is an SU(N) group element in the 
adjoint representation. The Weizsaker- Williams field is related to the valence color charge density 
by 



dib?(x) = f(x) 



(2.11) 



The classical field b has a well known diagrammatic interpretation. It sums the tree level diagrams 
to the one gluon component of the state created by the coherent operator C from the Fock vacuum 
(see Fig. 1). 




+ 





+ 



Figure 1: Diagrammatics of the classical field bf . The straight solid line represents the valence charges which 
serve as sources of the classical field. 



The explicit form of B has not been given in jy] , however the action of B on the gluon fields was 
calculated. Since we will use heavily this in the following, we give here the relevant expressions, even 
though they are quite lengthy. 



where 
and 



e(x) = - \9{x) - 9{-x)\ 



(Six) 



°° dp- 

„ 



d 2 q 



B(-x-)e*-* {t-l\ij{x,y)v^_ q {y) 



+0(x-)e l ^ x [T-L}%(x,y)v%_ jq ( y ) 



"P ,9 



+ h.c. , 



with 



Af b (x, y) = l Di—Dj + D l -^—D j - 2D i —d i 
13 y >yj \ l dD 3 Dd 3 3D J 



(x,y) ■ 



The covariant derivative is D bc = di5 bc — fuM and the projectors T, L, t, I are defined as 



L 



ab 
ij 



l D 2 3 



ab 



rpab raft jab. 

ij ~ "ij ij ' 



k 



d i^2 9 j' 



&ij lij ■ 



(2.12) 
(2.13) 



(2.14) 



(2.15) 



(2.16) 



The integral over p~ in eq. Q2.14 ) strictly speaking excludes the point p~ = since soft gluon modes all 
have non-vanishing light cone frequency. This subtlety however does not affect any of our calculations 
and we will therefore not indicate this explicitly. 

The transverse basis functions the analogs of plane waves and are normalized according 

to 

' d 2 qv;i q {x)vp q {y) = W£(x,y). (2.17) 
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where 



with 



(x,y) = f — V \^—[^5\x -y) + \c^{x,y)\ + -J—^S^x - y) - \c^{x,y)}\ 



1 



1 



C%{x,y) = \ 2d i —D j -2D i —d j 



ab 



(2.18) 



(2.19) 



For all practical purposes this means that for any finite frequency p the functions v can be taken 
as plane waves (we suppress rotational and color indices in which v can be taken as a unit matrix) 

1 



Vp-Jx) 



g iqx 



(2.20) 



For the mode at infinite frequency the normalization is different as is given by the operator Cfj{x). 
Although the infinite frequency mode does not contribute directly to any quantity which is pro- 
portional to rapidity, it does play an important role in insuring the completeness of the expansion 
basis in eq.( |2,12 ). We will not be working with this mode directly, but will simply use the fact that 
including this mode, the operator B is indeed unitary, with all the ensuing consequences. 

Since the color charge density is a quantum operator on the valence Hilbert space, it is also 
transformed by the action of B. This transformation has been found explicitly in flj. It turns out 
that for our present purposes it is not necessary to know the exact form of this transformation apart 
from the fact that it adds to j terms of relative order 0(a s ). 

We now consider the energy evolution of the hadronic wave function. Increasing energy of a 
hadron is equivalent to boosting it. Under the boost transformation, the longitudinal momenta of 



the soft gluons in the wave function eq^^) increase. As a result, soft gluons emerge from below the 
"cutoff" and the number of gluons in the wave function which contribute to any physical observable, 
such as a scattering amplitude increases (Fig. 2). 



Boost 



rapidity 
cutoff 



valence 
modes 



11 



t| + AY 



ay{ 



soft 
modes 



Figure 2: The light cone wave function before and after boost. Here n stands for gluon's rapidity, j for the 
color charge density of the valence gluons, and A for the soft gluon field. 



Thus in the leading logarithmic approximation the soft part of the wave function determines 
completely the evolution of physical observables with rapidity. However to find this evolution one 
has to make an additional step, and consider directly some physical observable. The purpose of 
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this paper is precisely to make this additional step from the soft gluon wave function to physical 
observables. 

The most basic and most inclusive observable is the forward scattering amplitude of our projectile 
hadron on some hadronic target. This is the observable we are going to discuss in this paper. It 
belongs to a broader class of observables which depend only on the correlators of the color charge 
density operator j a (x) and not on other characteristics of the projectile wave function. Evolution 
of observables of this type is governed by the functional evolution equation which generalizes the 
JIMWLK (21], |H and KLWMIJ |[§ - the Schroedinger equation of RFT. The kernel of the evolution 
equation can be viewed as the second quantized Hamiltonian of the QCD Reggeon Field Theory 
(RFT) as discussed at length in [|53| . By considering an arbitrarily dense hadron as a possible target 
in the scattering process, we will derive the functional evolution equation which is applicable in the 
full parameter range between the BFKL scattering of two small perturbative objects through to 
nucleus-nucleus scattering. The result of the present paper is therefore the RFT Hamiltonian which 
contains the effects of Pomeron loops. 

Another example of observables of this type are the single inclusive, double inclusive and in 
general multiple gluon production amplitudes, where all gluons are close to each other in rapidity. 



In the companion paper [56] we derive the explicit expressions for these observables in terms of the 
target and projectile fields, to be averaged over the corresponding wave functions, and discuss their 
evolution. 

We now briefly sketch what we have to do to derive Hrft- Consider the calculation of some 
hadronic observable O which depends only on the color charge density operator j a (x) and not on 
any other property of the wave function 



O = (v\ 6[j] \v) = I Dj W p [j] 0\j) . (2.21) 



The last expression is by now the standard way of representing the wavefunction average using 
functional integral. The weight functional W p provides a probability distribution for the valence 
charges j in the wavefunction. 

An important (but not the only) example of an observable discussed above is the eikonal S - 
matrix in the external field a 

S = expii / f{x)a a {x)\ . (2.22) 



In the approach of the expression of the scattering amplitude for a hadron-hadron scattering is 
given by 

S = ((v\S\v)) a (2.23) 
where the weight for the averaging over the fields a is determined by the wavefunction of the target. 

The same procedure applies to any observable. The matrix element O is the expectation value 
of the operator O in the projectile wavefunction. In order to obtain the final expression for the 
observable measured in the collision, we also have to average O over the target wavefunction. 

(O)t = [ DaW T [a] 6 (2.24) 
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where a represents all the target fields^ 



We assume, as before that only valence degrees of freedom contribute to the average eq.( |2.2lD 
at initial rapidity. When the system is boosted to rapidity Y, two important changes occur. First, 
the relevant wave function changes from \v) to \^f), and second the color charge density of the soft 
gluons has to be added to j a in the observable 0[j]. Both changes are due to the fact that the soft 
gluons after boost emerge with momenta above the cutoff rapidity, and now contribute to physical 
observables (see Fig. 2). Thus eq.( |2.21| ) at higher rapidity Y becomes 



Oy = <*| 0[f(x)+f soft {x)] |*> = (v\ &R a O[j] n \v) 
Here the charge density shift operator 



R„ 



fxisoft( X ) Sj a (x) 



and 



A a 

J soft 



(x) 



dk + 



a)\k + ,x) 



rpa I 

1 bc a 



■■(k+, 



(2.25) 



(2.26) 



(2.27) 



is the color charge density of the soft gluons with the integral over k + defined over the range of 
momenta corresponding to the boost parameter AY (Fig. 2). 

Eq.( 2.25| ) has a universal form, in the sense that it does not depend on the particular operator 
O. Since the state \v) is the vacuum of the soft gluon degrees of freedom, it is possible in principle 
to calculate the expectation value in eq.( 2.25| ) over the soft gluon part of the Hilbert space without 
the explicit knowledge of the operator O. This involves calculating the following matrix element 



(0| &[j, a, a f ] J- j °°f t{x) ™ £l\j, a, a 1 "] |0) 



(2.28) 



The quantum averages in eq.( p.28| ) are performed only over the soft gluon Hilbert space. Here the 
functional derivatives S/Sj are treated as c- numbers since according to eq.( [2.25 ) they should only act 
on j's in the operator O. The valence color charges j a in this expression are "almost" c-numbers, in 
the sense that they are not quantum operators on the Hilbert space of the soft gluons. One should 
however keep track of their ordering, since they do not commute with each other. 

The expression in eq.(2.28) of course depends on the total rapidity AY available to gluons in 
the soft Hilbert space. Since we are interested in deriving a differential equation in this rapidity, we 
have to treat this rapidity as small. To obtain the kernel of the high energy evolution equation we 
simply have to calculate eq.( 2.28| ) to first order in this total rapidity. 



lim 

AY^O 



(v\ ^R a O[j] O \v) - (v\0\j] 



AY 



= — Hrft 



H obi l«> 



where Hrft is defined as 
Hrft 



dY 



Y=0 ■ 



(2.29) 



(2.30) 



2 In most of our discussion below we will be focused on the evolution of the projectile and will not indicate target 



averages explicitly. One should remember, however, that the averaging over the target (2.24) is always assumed in final 
results. 
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We see that the evolution of any observable O is given by 

5 



±0 
dY 



(v\ H RFT 



3, TT 



0[j] 



Dj W y [j] Hrft 



3, tt 
°3. 



0[j] 



(2.31) 



Conforming with literature, in the last line we have represented averaging over the valence state as 
a functional integral over the color charge density with the weight functional W^fj]. The functional 
derivatives can be integrated by parts to make them act on W p . Assuming that the Hamiltonian 
Hrft is Hermitian (which we will find to be the case), we can rewrite this equation as an evolution 
equation for the weight functional as 



d 

dY 



Hrft 



Sj 



(2.32) 



Thus we conclude that in order to find the evolution of the scattering amplitude and other 
observables which depend only on the color charge density, we need to calculate the matrix element 
in eq.( |2.28| ) and expand it to first order in Y. This is what we will do in the next section. 

3. The Derivation 



We start by deriving a more explicit expression for the wave function 

|*) = n[j,a,a^} |0>. 



(3.1) 



We find it convenient to define the creation and annihilation operators in rapidity basis rather than in 
the basis of longitudinal momentum k + or frequency k~ as described in the introduction. Defining the 
rapidity variable as rj = In ^ with an arbitrary constant p$ , we rescale the creation and annihilation 
operators as 



(3.2) 




<{k + ,k) = ^ — a ?( V ,k) 

In terms of these operators we have 
"°° drjcPk 



Af(x , x) — -j= j 



2vr 4vr 2 



{f^ k \a,i,x-,x)a){r 1 ,k) + f< h ^ (3.3) 



with 



' Sij e 



- ikx 



(3.4) 



f (b t j,x-,x) = 

Strictly speaking the soft modes a{rf) live only at rapidities below Y, which is given by the parameter 
of the boost from \v) to |^). We do not denote this cutoff explicitly in most of our formulae and 
extend the rapidity integration over rj to infinity for all quantities for which the integration converges. 
We will indicate the Y dependence for divergent quantities and eventually this very dependence will 
determine the evolution in rapidity. 

The coherent state operator in rapidity basis is rather simple 



C 



cxp 



V2 



d 2 k 



drj 
27 



[aUri,k) +4 a (v,-k)} 



(3.5) 
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3.1 The vacuum of [3 



Given that 17 is a product as in eq.(2.7), we first analyze the state 

IO^sBIO). (3.6) 

Clearly, the state is annihilated by the operator (3 a = Ba a B^ . Here and in the following we 
denote all the indices of operators a, a' etc by a single Greek letter a. This includes rotational and 
color indices, as well as the transverse momentum and rapidity 3 . Since the operator B is of the 
Bogoliubov type, the transformation between the operators a, a) and (3, is linear homogeneous 
and can be written quite generally as 

p a = Ba a B ] = e a/3 a p + $ Q/3 oj , Pi = G*^ a\ + $* a/3 ap . (3.7) 

Since the transformation eq.( |3.7D is canonical, the transformation matrices Q and satisfy 

e$ T _ $e T = o, ee f - = 1. (3.8) 

The inverse transformation is given by 

«a = 9^ Pf) ~ *lp 0% 4 = Q T ap Pi - Qi p Pp. (3.9) 

This leads to another set of relations 

6 t G-$ T $* = l; G f $-$ T G* = 0. (3.10) 

The explicit form of the transformation matrices G and can be found from the results of [|J and 
we provide it later in eq.(3.37). For now however we keep the discussion general as the explicit form 
of Q and $ is not important until late in the game. 

To find the state \0)g we write it in the form 

lO)^ = F- l l 2 [K) e -| a - A -f4 |0). (3.11) 

Imposing the condition 

/3|0)/3 = (3.12) 

gives 

(6a + $at) e -| atAat |o) = (-GA + $)a t e-i aUat |0) = 0. (3.13) 

Thus we find 

A = Q- 1 ^. (3.14) 
To find the normalization of the state we have to calculate 

F[A] = (0| e -^ aAU e"5 at Aat |0) . (3.15) 



3 Below we will be frequently suppressing this index and using the more compact matrix notations instead. 
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For future use we calculate in Appendix A a slightly more general expression F[A, A] 



F[A, A] = (0\e 



■oA'a „ - i at Aa+ 



e - 3 atAaT = exp I — Trln(l - At A) 



(3.16) 



with an arbitrary matrix A. In the last expression the trace is over all the indices. Noting that 
F[A] = F[A, A] we finally find 



|0>/j 



e ^ 



|Trln(l-AtA) -|atAat 



|0) 



(3.17) 



with A given by eq. ( |3,14| ). Just like the classical field b, the matrix A has a diagrammatic represen- 
tation. It represents the two gluon component of the wave function obtained by action of B on the 
Fock vacuum, and can be represented by a two gluon vertex - Fig. 3. 






Figure 3: Diagrammatic representation of the two gluon vertex A. 



3.2 The matrix element. 

According to eq,( [2.30D we need to calculate the matrix element of the operator R a . Thus we need to 
multiply the state |0)/3 by the coherent operator C, eq.( [2.8D and then act on it with R a . The operator 
R a was introduced as the operator that shifts the valence color charge density by the color charge 
of the soft gluons. On the other hand when acting on the soft gluon operators it acts as the SU(N) 
rotation, rotating the soft gluon creation and annihilation operators by the unitary phase 

R\a b {x)R a = R bc (x)a c (x), R\a) b (x)R a = R bc {x)a) c {x) (3.18) 

with 



ab 



R ab {x) 



V exp{g I dt T 
Jo 



o Sj c (x,t) 



(3.19) 



As explained in to properly keep track of the ordering of the non-commuting operators j a (x) 
it is convenient to endow the valence color charge density j with the additional " ordering" coordinate 
t; j(x) — > j(x,t). For the purpose of our present derivation however, the only thing that matters is 
that R ab (x) is a unitary c - number matrix. 

Rotation of a and is equivalent to rotation of both indices of the matrix A as well as rotation 
of the classical field 6f (x) by the same matrix R(x). Thus we will consider the following general 
matrix element, which covers the particular case that interests us here 



Q[b - b] = (0| e~ 5 a « A L/3 a /3 e -ib a (a a +P a , 3 a J l3 ) e -ib a {a a + P a/3 aj J ) g - ^ aj, A a/3 ^ _ 

= (0| e -2 aAta e -«(ft-5)(a + Pat) e -i a tAat (3 2 q) 
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So far both the matrix A and field b are arbitrary. They will be specified shortly. We have introduced 
the matrix P a p to be able to reproduce eq.([3,5|) in the transverse momentum representation. The 
structure of the form f. b(k)[a(k) + a'(— k)] is reproduced in eq. fl3.20D by choosing the operator P to 
be diagonal in rapidity, color and rotational indices and to satisfy P(k,p) a(p) = a(—k). The operator 
P is clearly unitary and satisfies P 2 = 1. We will not need to know any additional properties of P 
in the following. 

The calculation is performed in the Appendix A with the result 

G[b -b}= F[A,A]exp|-i(6-6)(l-AP)[l-PA t AP]^ 1 (l-PA t )P(&-6)| (3.21) 
= F[A,A] exp|-i(6-6)(l-PA t )[l-AA t ]- 1 (l- AP)P(6-6)| . 

After these preliminaries we are almost ready to write down the answer for the matrix element 
in eq.( 2.3C| ). The only additional element we have to take into account is the ordering of the color 



charge density operators in the operators fit and fi. When calculating any physical observable, as in 
eq.( |2.25| ) all the color charge density operators in the operator O are ordered such that they are to 
the right of the operators j in fit but to the left of j's in fi. Thus when we introduce the ordering 
coordinate t, all j's in fit have to be assigned the value t = 0, while all j's in the operator fi must 
have value t = 1. In our derivation the fields b and A are associated with the operator fit and they 
depend on j(x, t = 0), while the fields b and A depend on j(x, t = 1). Also, as explained in |36|], since 
any physical weight functional W F depends only on powers of the matrix R(x) acting on 5[j], all the 
operators j that appear in the RFT Hamiltonian should be understood as acting on the matrix R as 
left or right rotations. 

f(x,t = 1) = gJUx) = "Str {^)T a ^^y| ; (3.22) 
f(x,t = 0) = gJ a L (x) = -gtr!.T a R(x) 



<5Pt( 



X J 

Jg(x) = [R(x)J R (x)] a / 
Thus we define two classical fields, bn and b^\ 

b a m = --f abc tf bd [J R ]d i XJ dc [J R ]= --r 6c 4 M ^C/f; 
b a Li = -- f abc U^ bd [J L ] d t U dc [J L ] =-- f abc ul bd di U dc (3.23) 
and similarly two matrices A 

A R = A[J R ] = Q-^Jr] <S>[Jr}; A l = A[J L ] = Q^[J L ] *[J L ] . (3.24) 



In fact for the purpose of the calculation of the matrix element of R a in eq. ( [2~30D the proper ordering 
of the factors of the charge density is equivalent to the substitution 

fit -». fit. =44; n^n R = c R B R (3.25) 
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where the subscript L (R) indicates that the respective operator depends on Ji (Jr). 

We are now ready to proceed with the computation of the matrix element fl2.28|) . We have 

(o| n[R a n R |o) = <o| n[n R |o), 



with 



= C L B L ; 



B L = RlB L , 



C L = R\C L R a 



(3.26) 



(3.27) 



The action of the operator Ql on the soft gluon vacuum is the same as the action of Q,l with the 
substitution 

b a Li ^b a Li = rf ab b b Li ; A L ^ A L = RjA L R. (3.28) 



Using eq. flA.10 ) and remembering the factor y/2 in eq. (|3.5| ) we can finally write 



(0| &R a n |0) = exp 
exp 

= exp 
exp 

= exp 
exp 

where we have defined 



~Trln(l - A* A) + ^Trln(l - A* A) + |Trln(l - A f A) j x 
(6 — &)(1 — AP)[l - PA^AP] _1 (1 - PA f )P(6 - &)} 

-^Trln(l - A 1 " A) + ^Trln(l - A f A) + ^Trln(l - A f A)| x 
-(& — &)(1 - PA f )[l - AAt] -1 ^ - AP)P(b - 6)} 

~Trln(l - A f A) + ^Trln(l - A f A) + |Trln(l - A f A)| x 
(b - 6)iV t [ee t - $$ t ]- 1 iVP(6 - 6)} 



(3.29) 



Nr = G r - <f> R P; 
Bl = &lR; 



$l = $lR; 



N L = G L -$ L P 
N L = N L R = @ L 



$lP. 



In eq.(3.29) we have omitted the subscript L on barred quantities and R on unbarred ones. This will 
be our practice also below whenever it cannot lead to confusion. 

There is one point which we have glossed over, namely the ordering of the factors of the charge 
density entering the classical field b relative to those entering A. It is obvious from our initial 
expressions that all factors of A in eq.( |3.29| ) have to be to the left of b, while all factors of A to the 
right of b; even though we did not indicate this explicitly in eq.( |3.29 ). As we will discuss shortly, 
however, this ordering is only important to subleading order in a s |TJ, and so we will not discuss it 
any further in the present paper. 

Eq.( |3.29| ) is far from being the end of the road. We now have to find A#, etc explicitly in terms 
of the matrix Ur (or equivalently in terms of the current Jr), and also expand the matrix element 
eq.( |3.29| ) to first order in the length of the rapidity interval Y. We will now address the first question. 
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3.3 The matrix A 



To calculate the matrix A consider expansion of the field Af(x ,x) in terms of two different sets of 
creation and annihilation operators a^,a and f3',(3. 



Referring to eqs. we can write expansion of A in terms of the operators 

Rfak) = Ba?( Vl k)BK (3.30) 

This has the form 

A?(x~ , x) = -L J°° g {gW"*) (a,i,x~ ,x) Pj(f],k) + g< b ^ (a,i,x~,x) (j* ( V ,k)} (3.31) 
with 



9r = 9 



^''"• fc )(o,i,x-,x) = (x\ Oi-x^e^* {l-2l) tj 8 ab + 



■El 

(x~) \e l 2k- 



1 - 2L R )]$ + e(x") [A fl (l - J - (3.32) 



where as usual \k) is the eigenstate of transverse momentum and \x) is the eigenstate of transverse 
coordinate. The subscript R in the above expressions indicates that the relevant quantity depends 



on J R . Here L$ R = Df{b R 



ed 



Df(b R ). 



The two sets of functions (/, /*) of eq. fl3.4|) and {g,g*) of eq.( |3.32|) constitute complete bases on 
the functional space. It is convenient to define a scalar product on this space as 



(g\f) 



[(c,n,£,p);(b,j,r),k)] 



*| f\[( C > ri ^:P);( b J:'?. fc )] 



(5*1/) 



J dx-d 2 xd + g*^ n ^ p \a,i,x-,x)f^' j ^ k \a,i,x-,x); 

a,i 

J ' dx-d 2 xd + g^*\a,i,x-,x)f( b >>> r i> k \a,i,x-,x). (3.33) 



With this scalar product the bases defined in eqs.(3,4. 3,32j ) are orthonormal: 
^(b,j,r,,k)\f(eM,P)) = s^SjiSiri-^ik-p); ( g Q>j,i>Q\gW£ffi)) = SPSflSb-OPik-p) . (3.34) 

The resolution of identity in this space (completeness relation) is given by 
dr\ d 2 k 



2vr 4tt 2 



f( a ^ k )(b,j,x-,x) d+ f< a ^ k \c,l,y-,y) 



(3.35) 



f< a ^ k \b,j,x-,x) d + f^ k \c,l,y-,y)] = 5 bc 6 jl 5(x- -y-)6 2 (x-y) 



similarly for the basis (g,g*). 

We can now calculate the matrices and $ defined in eq.(3.7) 



6 = (g\f); 



* = ( 9 \n 



(3.36) 
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Performing the integration we find explicitly (see Appendix B) 



9(p, rj; k, = i(p\2(l -L R ) + (l-l- L R )A R \k) + 

: (p|(l-2Z)|A;>-i(p|(l-2Lfl) ' 



U R 



Hp, m -k, = i(p\2(i -l r ) + (i-i- L R )AUk) + 



^v(p, »7; fc, = ©(p, »7; k, - $(p, ??; -fc, 



■<p|(l-2Z)|A ! )-(p|(l-2Lj l ). 2 ' 



(3.37) 

In the last terms in all three expression the ordering of the factors is such that all d 2 are to the 
right of all D R . We have also not indicated explicitly the color and rotational indices on G and 
Since the basis function g in ( [3.32 ) depends on J R , all the function B, $ and N computed above are 
"Right" quantities and should be understood as Q R , & R and N R . Substituting Jl for J R generates 
the analogous "Left" quantities. 

Mindful of the derivations in the previous subsection we find it convenient to define yet another 
set of creation and annihilation operator (3 and /3' by 

fi a = Ba a & = Qvpap + $ aP a\, 4 = ©a/3 a P + K/3 4 • ( 3 ' 38 ) 

Using eq.( |3.7D we find that the set /?, is related to (3, ft by 

p a = (eet - $&) af) p p + ($© T - e$ T ) Q/3 % (3.39) 

We thus see that the basic matrix that enters the last line of eq.( [3.29 ) K = 00^ — is the 
transformation matrix between the operators (3 a and /3 a . We can obtain it by calculating the overlap 
matrix between the appropriate basis functions 

K a(3 = (00t - $#t )a/3 = ( 5a |^) 

£ Q/3 = ($© T - ©$ T ) Q/3 = (g a \9p) ■ (3.40) 

Although we do not have the explicit form of g, we can still use eq.( |3.40| ) to calculate the matrices K 
and E using the following argument. We do have the explicit form of the basis functions gi, which 
are the same as g R except for the substitution J R — » J^. We can thus calculate explicitly the overlap 
matrices between g R and gz, (the analogs of K and E). We then notice that the operator (3 is related 
to the operator (3l by (3 = R^(3l, and thus ©l = @lR', &l = &lR- Then referring to eqs. (|3.40|) 
we see that to calculate K and E from the overlap matrices of g R and gi we must rotate all terms 
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involving 6^ by the operator w on the left. Performing this calculation explicitly we find 



K(p, rj; k, £) = i(p|2(l - Z - L fl )ljt(| — L£) + (! — / — L R )(A - &l)R\l -l-L L )\k) 



+ 



1 - e*-" + 



: <p|(l - 2Z)i?t(l - 2/)|fe) - i{p|(l - 2L fl ) 



1 - ^ef-" - ze 



flt(l-2£ L )|fc) 



S(p, 77; —k, = *<p|2(l - Z - LiO# f (Z - L x ) + (1 - I - L R )(A - A R )R\l - I - L L )\k) 



+ 



1 + e?" 7 ? 



— (p|(l - 2Z)i?t(l - 2/)|fc) - i(p|(l - 2L R ) 



1 + ^fe^ 



— 2ix,)|fc) 



with 

D = R^D L R; A = R ] A L R. 

In eq. ( |3.4l|) all Dr are ordered to the left of all Dl. 

Note that similarly to G and the matrices K and E satisfy 



(3.41) 



(3.42) 



KE T — EK T = 0, KK^-EE^ = 1; K^K — E T E* = 1; K^E-E T K* = 0. (3.43) 

For future convenience we also define the analog of the matrix A: 

S = IT^S. (3.44) 
Using the identities eq. (gg jjTjjlgjg ) it is straightforward to show that 

= J&K ; = K T K* . (3.45) 



4. Hrft 

Now that we have all the ingredients, we are ready to put them all together and to derive the RFT 
Hamiltonian. Before continuing with this calculation we want to clarify the counting of powers of a s . 

4.1 Counting powers of a s 

First, using eq.( 3.45| ) we have 



•~Trhi(l-AU) + ~Trm(l-AU) + -Trln(l-A t A) = ~Trln(l- 



--TrlnfMt). (4.1) 



Now, recalling the definition (2.30) and referring to eq.(3.29|) and eq.(|4.1|) we can write 



H 



RFT 



dY 



(b-b)N* K~ l N P(b-b) + - TrlnKK^ 



y=o ■ 



(4.2) 



The two terms in this expression are not the same order in a s . We remind the reader that this 
expression was derived starting with the evolution of the wave function of p]]. This evolution kept 
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the leading in a s term for arbitrary parametric value of j and also terms which can be 0(1) at large 
j. The former terms in the evolution of the wave function give rise to the first term in Hrft eq.( |4.2| ), 
while the latter leads to the last term in eq.( |4,2D . This can be seen directly from eq.(4.2). As long 



as j <C 1/g, the classical fields 6 and 6, as well as the operators A and A can be expanded in powers 
of gj. To count powers of a s we have to decide how to treat R. We know that when we calculate 
the forward scattering amplitude, after averaging over the projectile wave function, R turns into S - 
the S-matrix of scattering of a single projectile gluon on the target [ 31 ],E3]. Therefore if the target 



is dense we should count R as being of order one but not parametrically close to the unit matrix. In 
this case we see that 

(b - W ~ 0(j 2 ); KK^ = 1 + 0(g 2 j 2 ) . (4.3) 

Thus indeed the first term in eq.( |4,2| ) is leading, while the second is suppressed by a s . 

When j ~ 0(1/(7) the matrix A can not be expanded. Still A ~ O(l) while b — 6 ~ 0(l/a s ), and 
the first term is leading. 

If the target is not dense, the matrix R is close to unity R = 1 — 5R. Then b is perturbatively 
close to 6; that is (6 - b) ~ 5Rb, and also K ~ 1 + 0(5R). The first term in (|]|) is (5Rb) 2 , and 
dominates over the second term, which is (5R) 2 . 

Thus we conclude that in all parametric regimes of the projectile and target fields, the second 
term in eq.( |4.2j ) is suppressed by a s relative to the first term. In most of this paper we will therefore 
neglect the last term and will only consider the leading order piece 



H RFT = ± 



(6-6) N ] K^NPip-b) 



\Y=0 ■ 



(4.4) 



At this point we are ready to comment to the question of ordering of the factors of 6 and A briefly 
mentioned in the subsection 3.2. Once we restrict ourselves to eq.(4.4), this ordering becomes irrel- 



evant. As we have noted above, in principle all factors of A have to be ordered to the left of 6, etc. 
However, changing the order of A and 6 brings a term proportional to the commutator 



The additional terms this procedure generates are of the same order as the last term in eq.( [4,2[ ) which 
we have neglected. Thus the ordering question is irrelevant as long as we are interested in the leading 
order term. 

The last term in eq.(|4.2j) originates directly from the Bogoliubov part of the transformation 
S7. However, even though this term itself in eq.(4.2) is suppressed by ct s , it is important to realize 



that the Bogoliubov operator contributes in an essential way to Hrft and to other observables at 
high energy. For example the factors of N, K~ l and N in eq.([4.4|) are direct manifestations of the 
Bogoliubov part of the transformation Vt. The presence of the Bogoliubov operator B also leads 
to interesting effects in other observables. For example, as shown in the companion paper |j56| , it 
is responsible for the appearance of short range rapidity correlations in multi-gluon spectrum. For 
these particular correlations B gives the leading effect, as these correlations would be absent if the 
transformation $7 had only the classical piece of the type of C. We discuss this in more depth in 
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4.2 JIMWLK for the umpteenth time 



At this point to keep our finger on the pulse it is a good idea to perform a consistency check of our 
calculation so far by reproducing the JIMWLK equation. The JIMWLK equation is obtained by 



expanding eq. ( |4.4j) to second order in 5/dj. Consider the expression 

(b-b)N* K~ l N P(b-b). (4.6) 

Since b — b is already first order in the derivative, in this approximation we have to take A = A and 
N = N. This simplifies things dramatically, since in this case K = 1. Thus in this limit we only 
need to calculate 

ib-l) / N\ n ,i)N{U)Pib-l) (4.7) 

where for simplicity we do not indicate the transverse coordinate dependences. Note that here and 
in the following the measure of the rapidity integration is defined as drj/2-K, so that J = J 

Since N(rj,^) = N(r) — £) (see ( |3.37D ) is the function of rapidity difference only, one of the rapidity 
integrals in eq.( ^4.7[) gives the factor of Y, and the RFT Hamiltonian becomes 



Hrft = ^-(b-b)N[N ± P(b- b) (4.8) 

with 

N ± = I N(rj) . (4.9) 



n 



We remind the reader that b is independent of rapidity. The integral over the rapidity in both terms 
of N in eq. ( |3^37|) is the same (after the appropriate shift of the integration variable) 

/ dr] - = P I dri +2ivr / di]5(e v - e~ v ) = in. (4.10) 

J e v — e V — ie [J e r] — e ri \ J 

Thus we obtain 

iVj_ = [1 — l — L]. (4.11) 

With this result, expanding b — b to first order in 5/5 j in eq.( |4.8j ) it is straightforward to see that 
eq.( f4.8p reproduces the JIMWLK Hamiltonian in the form obtained in M. Having established this, 
we now turn to the calculation of the Holy Grail - the complete Hrft- 

4.3 The Holy Grail: the RFT Hamiltonian. 

We have to evaluate 

[b - b] R- 1 NP [b-b] = 

[b(x) - b(x)} [ N^ V ,x;Z,y)K- 1 (Z,y;\,z)N(\,z;(,u) [b(u) -b(u)] (4.12) 

x,y,z,u Jr],t;,\£ 

Note that all the operators in question (N, N and K) depend only on rapidity differences. Thus the 
integration over rapidities yields 

1 



2tt 



[b(x) - b{x)\ N[(x, y) Kl\y, z) N ± (z, u) [&(«) - b(u)\ Y . (4.13) 

x,y,z,u 
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Here Y is the total rapidity interval and we have defined 

/dfj 
— K(t],z\£,u) 



(4.14) 



Note that since the operator K depends only on rapidity difference, the integral of the inverse is the 
same as the inverse of the integral 



(4.15) 



If all rapidity integrals were finite, the expression in eq.( 4.13| ) would be proportional to the total length 
of rapidity interval opened by the boost, or "created" by the action of the operator f2. Eq.( |2.30| ) 
then would imply 



H RFT = - 



[b(x) - b(x)] NUx, y) K~ L (y, z) N ± (z, u) [b(u) - b(u)] 



i x,y,z,u 

Using eqs.( [4.1C' , 4,ll| ) we have 



N ± (x,y) = l-l-L R ; 



N±(x,y) = (l-l-L L )R. 



(4.16) 



(4.17) 



The integration of K 1 however gives a puzzling result. If expression eq.(3.41) is taken literally, we 
find that at large Y 

[ K(ti,£)<xY; I K' 1 ^) cx 1 (4.18) 

It then seems that eq.( |4.13j ) does not have the overall factor Y and thus does not yield a logarithmic 
evolution. The resolution of this puzzle is the following. Since Hrft is obtained by differentiation 
with respect to Y, at this stage of the calculation we really should consider Y to be small rather than 



large. In writing eq.( 4.13 ) we have assumed that it does not contain higher powers of Y. On the other 
hand our expressions when expanded in powers of a s at small j contain higher order perturbative 
contributions. The standard perturbative expansion at any fixed order (beyond the leading one) 
contains higher powers of Y. The simplest source of higher powers of rapidity is the "iteration" of 
the leading order kernel. Such terms have been identified and removed in the NLO computation 
of the JIMWLK kernel ^7|. Obviously in our case also these higher powers of logarithms have to 
be subtracted in order to extract the evolution kernel. This of course is the same as treating Y as 
infinitesimal while extracting Hrft from the fixed order perturbative calculation. Clearly in the 
present case we should do the same - we should learn how to subtract the higher powers of Y in 
eq. (fDJ). That such higher powers exist is quite obvious. If we think of Y as of finite range of 
integration of the rapidity integrals, eq,( 4.18| ) reads schematically as 



K(v,0 = «l + k 2 Y; 



,. CO 

Jr > n=0 



(4.19) 



with K\ and k 2 some operators in the transverse space. This structure is demonstrated explicitly in 
the Appendix C, where we perform the perturbative expansion of eq.( 4.13| ). We show there that the 
next to leading term contains a term proportional to Y 2 , as in [47]. 
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In view of this it is clear that our task is to extract the n = term in the sum in eq.( |4.19| ). 
This is actually not difficult. Examining the expression eq.(3.41) it becomes clear that the K2 Y term 
in eq.(4.19) comes from the integral over rapidity of the rapidity independent terms in K. Let us 
rewrite K as follows 

K(ri - = K{v ~ + K A {rj - £) + K c (4.20) 

where if is a symmetric (K(r] — £) = K(£ — 77)), bounded operator whose matrix elements vanish as 
\r] — £| — ► 00; while Ka is antisymmetric {Ka(j] — £) = —Ka(£ — ??)) and Kc does not depend on 
rapidity. We can then write 



/ k( v -0 = [ K(n-0 + k c y. 



(4.21) 



Thus expansion of J^K^ 1 in powers of Y is the same as expansion in powers of Kc- Since we are 
interested only in the leading term in this expansion, we need to drop Ka + Kc from the expression 
for K in eq.(3.41) and retain K. Eq.(OI) now reads 



H RFT = - 
A quick calculation gives 



[b(x) - b(x)] Nl (x, y) Kl l (y, z) N± (z, u) [b(u) - b(u)} . 



(4.22) 



x,y,z,u 



K(V,0 = i(l-2Z)i?t(l-20 



2(1- 2 ^) 



1 



1 

To calculate the integral over rapidity, consider an integral 
1 1 



+ 



+ 



1 



IC 



1 



D 2 



R\l - 2L L ) (4.23) 



dx 



1 



+ 



it 



1 



1 



ae 



IE 



dx- 



1 



ae' 



ie)(l 



ae~ 



le 



(4.24) 



To integrate we can close the contour in the complex plane. The integrand decays exponentially 
at infinity everywhere except close to the imaginary axis. Nevertheless the contour can be closed. 
To see this note that the integrand has an infinite number of poles at x~ = — In a + 2irin — ie and 
x£ = In a + 2irin + ie. The residue of the pole x~ is the negative of the residue of the pole x+. Thus 
the contribution of a pair of poles vanishes, and the result of the integral does not depend on how 
many pairs of poles the contour encloses. Closing the contour above the real axis we see that the 
only uncanceled pole is Xq , it's residue is unity, and the result of the integration is 

1 



Analogously 



dx 



dx 



1 



1 



ae- 1 



if 



+ 



1 



1 



+ 



ae 



1 



ie 



ae- 1 



ie 



1 — ae x + ie 



2iri . 
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Putting it all together we find 

1 

2 



K± 



(1 - 21) (1 - 21) + (1 - 2L R ) R ] (1 - 2L L ) 



(4.25) 



(4.26) 



(4.27) 
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Thus the complete RFT Hamiltonian is given by 



Hrft = -[b R tf-b L ] (l-l-L L ) 



TT 



(1 - 21) i?t (1 - 21) + (1 - 2L R ) ^ (i _ 2L L ) (I- I- L R ) [b R - b L ] (4.28) 



This is the main result of the present paper. Note, that Hrft can be written entirely in terms of 
three unitary matrices R(x), Ur{x) and Ul(x). Restoring all the indices and transverse coordinate 
dependences we have 



Hrft = ^3 / [b Ri {x)R^(x) - b%{x)] 

Jx,y,z,z 

x ~5 ac + [Ul(x)U L (z)] ac ] K~j k cd (z,z) 



1 



13 (x — z) 2 

1 



(x - z)j (x - Z)j 
(x - z) 4 

, (y - z) k (y - z)i 



(y - z) 2 

5 de + [uUz)U R (y)] de ] [b%(y) - R^(y)b[ k (y)) 



(y ~ z) 4 



(4.29) 



with 



2vr 2 J z 



(x - z)< 



(x - z)j (x - z) h 
(x - zY 



(z - y) k (z - y)j 



{z-yf 



{z-yY 



x \rW(z) + \uUx)U R {z)R\z)uUz)U L (y) 



ab 



and 



°L(R)i 



.1 fabc u W djU dc 
9 



'L(R) 



L(R) ■ 



(4.30) 



(4.31) 



We note that in appropriate limits eq. fl4.28l ) reduces to the known results. We have already 
discussed the JIMWLK limit, and it also directly follows from the above formula. The leading order 
expansion of R in powers of S/Sj requires that we set R = 1 everywhere except in biR, which has 
to be expanded to first order. In this limit eq.( [4.28 ) reduces to Hjimwlk- In the limit of the low 
charged density on the other hand, we have to expand in j. To lowest order in j we have Lr = Ll = I. 



The denominator in eq.( 4.28| ) cancels the two adjacent factors (1 — 21) in the numerator leaving the 
factor 1/2. Additionally we have to expand b to first order in j, that is b{ = We thus obtain the 
KLWMIJ limit lETO 



Hrft — > Hklwmij = — 

ZTT 



3%-jR^R 



d 2 



d 2 



d 2j d 2 3 



(4.32) 



Note that if we do not expand b to first order on j but still set Lr = Ll = I in the rest of the 



expression, we obtain the form dubbed in [31] "KLWMIJ+" 



Hklwmij+ = 7; [bR - R^bi} 2 

ZTT 



(4.33) 



Although this looks simple and rather appealing, it is not a leading term in Hrft m any well defined 
limit. This expression was also derived in [^, 42 j. It has been suggested in |5l] that it might in fact 
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be the complete expression for Hrft-, including the Pomeron loop effects. Our present derivation 
makes it clear that this is not the case, as in general the terms kept in Hklwmij+ ar e as important 
as the terms omitted. In particular none of the effects coming from the Bogoliubov part of the 



operator Q is taken into account in eq.(4.33) and thus the JIMWLK limit is not reproduced . 
5. Discussion. 

First, let us summarize the main results of this paper. We have calculated the Hamiltonian of the 
Reggeon field theory in the eikonal approximation. 

Hrft = -\b R - b L R] i?t (1 - I - L L ) (5.1) 

7T 

;i - 21) R) (1 - 21) + (1 - 2L R ) i?t (1 - 2L L j\ ~* (1 - / - L R ) [b R - b L ) . 

This Hamiltonian governs the evolution of hadronic observables with energy. It reproduces the 
JIMWLK and KLWMIJ limits discussed in the literature earlier, as well as the "intermediate" form 
- KLWMIJ+. The evolution generated by this Hamiltonian is valid for all interesting values of the 
color charge density g < j < ~. This covers the hadronic targets from a dilute perturbative object 
(a "dipole") to a dense non-perturbative object with large gluon density (a "nucleus"). Importantly 
it also describes the evolution of a nuclear projectile at all momentum scales. This is important 
since at high transverse momentum - larger than the saturation momentum Q s , a nucleus is also 
transparent, and thus the color charge density at these momenta is small. Thus if one is indeed 
interested to follow evolution at large momentum scales, the JIMWLK equation is not adequate even 
if the projectile is a nucleus. Another regime where JIMWLK is not adequate for a nuclear target 
is in description of nuclear periphery where the color charge density is also small. Again, in this 
regime Hrft supercedes Hjjmwlk, although here the situation is complicated by non-perturbative 
soft effects related to confinement physics. 



We note that the Hamiltonian eq. (|5,l| ) is invariant under the Gribov's signature symmetry, which 
in [53] was identified as the symmetry under the transformation R —* R* , Jl —> Jr- 



An attempt to derive the complete Reggeon Field Theory Hamiltonian was made in [ 38 1 . As 



we mentioned earlier, the result presented in [38] is the KLWMIJ+ Hamiltonian, while the correc- 



tions due to the Bogoliubov part of the operator Q were missing. Since the derivation of [38] is in 
the Lagrangian formalism, it also leaves unanswered the question about the canonical structure of 
the operators entering the RFT Hamiltonian. Our derivation clarifies this question along the lines 



discussed in [54]. The Hilbert space of the Reggeon field theory is the space of the functionals of 
the unitary matrix R: ty\R{xj\. The unitary matrix R is therefore the only independent quantum 
degree of freedom in this theory. The operators Jl{%) and Jr(x) act on this space as the generators 



Ref.[|38j asserts that eq.(4.33) does reproduce the JIMWLK limit. The argument of |3q ] however is not direct but 
is rather based on an application of the Dense-Dilute Duality transformation J35| Ur — > R. To be able to apply this 
transformation to Hrft however one has to know the transformation properties of Ul ■ A certain transformation of Ul 
was postulated in |38). Unfortunately this did not take into account the fact that Ul is not independent of Ur and R, 
but rather is unambiguously determined once Ur and R are known. The transformation postulated in turns out to 
be inconsistent with this dependence, and thus the argument about reproducing the JIMWLK limit given in feq] fails. 
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of local left and right SU(N) transformation group. Consequently the operators Ul and Ur do not 
commute with R, but have rather complicated commutation relations determined by the dependence 
of Ur on Jr and of Ul on Jl- On the other hand Ur commutes with Ul, since Jr and Jl commute 
with each other. 

There is a question, which we have not touched upon so far, but which is important for under- 
standing the consistency of the whole approach. The issue is when can we consistently treat the 
length of the rapidity interval Y as large, and when can we treat it as small. In our derivation we do 
both. First, when deriving the evolution of the hadronic wave function in pj, we have approximated 
the interaction between the soft modes and the valence modes of the gluon field by the eikonal ver- 
tex 5 . This is valid as long as the main contribution to physical quantities comes from the soft gluons 
with rapidities much smaller than the rapidities of the valence ones. The borderline region, where the 
rapidities of the "soft" and "valence" gluons are comparable, is of order 5Y ~ 1. Thus consistency 
of the eikonal approximation requires F>1. As long as the eikonal approximation is valid, we can 
ignore finiteness of Y and treat the phase space available to soft gluons as infinite, which is what was 
done in Q. 

On the other hand when calculating Hrft itself we have differentiated the matrix element with 
respect to Y at Y = 0. Thus at this point we have treated Y as small. Even though it may seem 
odd, in principle the two approximations are not necessarily incompatible. The question is how fast 
is the evolution of a given physical quantity. For example, in the BFKL calculation to get significant 
change in the cross section due to the contribution of soft gluons, one needs to increase the rapidity of 
a process by Y ~ This means, that when taking the derivative to derive the BFKL Hamiltonian 
the smallness of the step in Y means simply Ay <C — . At small coupling constant this is still 
compatible with the condition for validity of the eikonal approximation Ay 3> 1. We can thus use 
the eikonal approximation ("large" Y approximation) to calculate the evolution even on a "small" 
rapidity interval. One can have one's cake and eat it! 

The same exact argument holds for the KLWMIJ and JIMWLK Hamiltonians. Formally the 
evolution eigenvalues in both cases are of order a s . To see this remember that rescaling j, the 
KLWMIJ Hamiltonian can be written as 

Hrlwmij = ct s Hklwmij[~3, R] ■ (5-2) 

9 

We use this Hamiltonian when j ~ g and R ~ 1, thus parametrically the eigenvalues are ujk lw mij ~ 
0(a s ). Likewise for JIMWLK 

H JIMWLK = Oi s HjimWLk[U, —, — rr] (5.3) 

where U ~ 1 and gj ~ 1 so that u jimwlk ~ 0(a s ) 6 . Thus for the evolution of the forward scattering 
amplitude in either KLWMIJ or JIMWLK limits, we have 

S(Y) ~ e asKY 5(0) (5.4) 

J This is technically not the same as the eikonal approximation for scattering on the target, although the consistency 
of the whole approach likely requires both. 

6 In fact we know that due to the dense-dilute symmetry p3 the two sets of eigenvalues are equal. 
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where k is of order one. Thus again we see that the S'-matrix changes significantly only if the 
rapidity is large Y ~ ±. So in the KLWMIJ and JIMWLK limits one can have large rapidity 
interval to accommodate the eikonal approximation, but still small enough steps in rapidity to be 
able to calculate the derivative. 



This argument cannot be straightforwardly extended for Hrft of eq. fl5.ip . Assuming b — b ~ 
0(1/ g) and R ~ 1 we expect Hrft ~ 0(l/a s ), and thus lorft ~ This suggests that the evolution 
with Hrft is very fast and observables change significantly over the rapidity range Y ~ 0(a s ). If this 
were indeed the case, the whole approach beyond the KLWMIJ/ JIMWLK limits would be in question. 
However, the situation is not really so bad. Although we do not know the spectrum of Hrft it is very 
likely that just like in the JIMWLK/KLWMI J case the spectrum is continuous with lowest eigenvalue 
being zero. In fact the existence of eigenstates with zero eigenvalue is straightforward to establish. 
Just like for JIMWLK/KLWMIJ the state (in the RFT Hilbert space) with zero charge density 
j|Yang) = and the state with wave function independent of j: i?|Yin) = |Yin) are both eigenstates 
of Hrft with zero eigenvalues |54|]. Physically those are the totally white (vacuum) and the totally 
black (black disk) states. If indeed the spectrum of Hrft is continuous then the spectrum also has 
eigenvalues of 0(a s ). Evolution of observables which have significant overlap with eigenstates that 
correspond to these eigenvalues then has the same speed as for the JIMWLK/KLWMIJ evolution. 
More formally we can represent any observable and its evolution as 



(O) = j[dj]W p \j] 0[j] = (W p \0) (5.5) 
with 0[j] is itself a target average of the original observable operator 

0[j}= J[dS]0[j,S]W T [S}. (5.6) 
and bra and ket referring to states in the Hilbert space of RFT. Then 
^(O) = - I [dj}W p [j]H RFT 0[j] = -Y,"n(W p \n)(n\0)- (0) Y = £ e~^ Y {W p \n)(n\0) . 

J n n 

(5.7) 

Thus if at the initial rapidity both \W P ) and \0) have overlap with eigenstates n such that u! n ~ 



0(a s ), the evolution of such observables with Hrft is slow. In the companion paper 56] we discuss 
an argument purporting to establish that inclusive multi-gluon amplitudes evolve with the JIMWLK 
Hamiltonian. If true, this would then suggest that the multi-gluon amplitudes are indeed observables 
of this sort. 

An interesting question is whether the forward scattering amplitude in the nucleus-nucleus col- 
lision is also such an observable. As we have discussed above, the forward amplitude is not evolved 
with JIMWLK Hamiltonian, but with complete Hrft- Also if we consider scattering on a fixed 
configuration of the target field a, the evolution is certainly fast. The counting which gives Hrft as 
0(1/ a s ) clearly holds for a generic configuration of a ~ 1/g. A small change of the projectile charge 
density 5j ~ 0(1) leads to a large change of the phase since a is very large. Another way of putting 
it, is that the observable exp{ija} is orthogonal to |Yin) and probably also almost orthogonal to 
other Yin -like states, while being the weight function of a nucleus, has a very small overlap 
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with |Yang) and Yang -like states. Thus the overlap of W p with expjij'a} is dominated by the states 
n with large eigenvalues u n . 

The interesting quantity is however not the scattering amplitude on a fixed configuration of the 
field a, but rather the target-averaged observable like in eq.( |5.6[ ). Averaging over the target fields 
has the effect of averaging the phase expjij'a} to an extremely small value. Thus it could well be 
that the contribution of the states with large uj n is erased by the averaging procedure. If this is 
the case we are back to the situation where the states with small u n dominate both the value of 
the observable and its evolution. In this case the eikonally derived Hrft is applicable also for the 
forward scattering amplitude. We have not studied this question any further and at the moment 
cannot make a definitive statement. 

Finally we want to mention several possible lines for further research. There are two questions 



that can be addressed immediately. First is extending our calculation [56] of multi-gluon spectrum 



to the situation when the rapidities of the observed gluons are far from each other in rapidity. This 



should follow very closely the derivation of [51| and should not pose any major complications. 

Secondly, we have not studied here the effects of the Trln term in Hrft, since it is a sub leading 
in the coupling constant effect. At j < 1/g this term collects some of the subleading corrections, but 
may not be complete (we have not analyzed this question). However, as follows from jjj], at j ~ 1/g 
this term contains all the NLO corrections to Hjimwlk- Thus for large projectile fields it is very 
interesting to calculate this term explicitly, as it is the generalization of purely gluonic terms in the 
calculation of 47] for the case when the projectile is dense (and the target, in principle is arbitrary). 



A more formal question is that about the selfduality of Hrft- It was shown in [^] that the 
RFT Hamiltonian must be selfdual under the Dense Dilute Duality transformation Ur <-> R. We 
have not addressed this question in the present paper. The selfduality is not an entirely trivial issue, 
since to verify it we need to understand the duality transformation properties of Ul- We are hopeful 
that this can be done and plan to address this question in the future. 

Finally, of course it is imperative to understand what physical consequences has the inclusion of 
Pomeron loops in the evolution. We hope that it will be possible to analyze the dynamics generated 
by Huft a t least numerically. It appears to be a much more complicated problem than, for example 



JIMWLK evolution successfully studied in the Langevin formulation in [57], as the kernel now has 
an infinite number of derivatives and does not allow Langevin formulation. Still we are hopeful that 
some quantitative analysis will be possible. It is generally true, that it is one matter to find the Holy 
Grail, and quite another to learn to drink from it. 
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A. Appendix - The calculation of matrix elements. 

In this appendix we give the details of the calculations in Section 3. We have to calculate the matrix 
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element 



F[A, A] = (0\e--2 aAta e -^ atA ^\0} 



(A.l) 



Differentiating this with respect to A we obtain 



dF 

dA a/3 



(0|e- 



r (iAa 



1 f f 



e -! atAat |0) 



= ( |-^-^- e -^ Aa e-| atAat |0) 



9a Q Sag 



<0|e" 



a A a 



2 a/3 

The resulting differential equation: 



1 - r) F - 



dF 
9At J 



9A, 



a/3 



9At' 



is solved by 



F[A, A] = exp<j --Trln(l 



At A)) 



(A.2) 



(A.3) 



(A.4) 



We now turn to the matrix element 



G[ la ] = (o\ e -h a ^ a e - fb( - a+Pat h- ib ( a+Pa ^e-^ aiAat \0) = (0|e _ 3 
where we have introduced notation 

7a = b a b a . 

To calculate this matrix element we differentiate with respect to j a . 
dG 



\ a A* a -17 (o+Pat) at Aa+ 



e 2' 



|0) 

(A.5) 
(A.6) 



^7a 



i(0|e-2 aAt a (o a + P^aJ) e -^ a + Pa+ ) e"l fflt Aat |0) 



i(0| 



a/3 



e 2 J 



-i7(a+Pat) e -|Aot a 



|0> 



i (S a p - (PAt) Q/3 ) (0|e-5 ° At a 0/3 e-^( a+Pat ) e -5 at Aat |0) 

<W - (P^U) k7P)^+(0|e-^ At - e "n( a W) ( _ A/37a t )e -iatAat |0) j (A ?) 



The last term can be expressed in terms of J^- if in the very first line instead of commuting to 
the left, we commute a to the right 



dG 



i(0\e 



■oAt« 



+ Papal) e-*T(°+ Pot )e-3 at Aat |0) 



= i(5 Q/3 -(AP) a/3 )(0| e ^^ Ata e-^ a+Pat )(Pat) /3e -^ tAat |0) 



So that 



<0|e-2 aAt a e- i ^ a+Pa ^ ) ai e"3 fflt Afflt |0) 



P- 



1 - AP 



a/3 97,3 ' 



(A.8) 
(A.9) 
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Together with eq.( |A.7l ) this gives equation 



_ = _(!_ PA^apPwuG - (1 - PP) a p[AP]p e (l - AP)7J^- 



This is solved by 



(A.10) 



(A.ll) 



G[i\ = G[0] exp{-^ 7 (l - AP)[1 - PAUpJ-^ 1 " P^)Pl} 
= G[0] exp{-^ 7 (l - PA f )[l - AA + ] _1 (1 - AP)P 7 } 

where G[0] = F[A, A] of eq.(|3l^). 

B. Appendix - The overlap matrices. 

In this Appendix we calculate the overlap matrices G and <3? used in the text. Here we will use 
frequency p~ rather than rapidity to label the basis functions. However the normalization is still 
taken to be consistent with expanding the field A in terms of creation and annihilation operators 
a„, a\. We start from the expression for the basis functions 



(B.l) 



[a,i,x ,x) 



(x\e(-x-)e^ x (1 - 2l) ij 5 ab + 9(x-)[e l ™= x (1 - 2Z,)]# + e(aT)[A(l — I — L)\ 



Thus 



8+g*P = 5(x-)(p\2l - 2L + (1 - I - L)^\x 



ab 



(B.2) 



2p~ 



(p\9(- X -)e %>- d\l - 21) + 0(aT)[e D\l-2L))\x) 



In the expression above we have not written all the indices explicitly. We first check the orthonor- 
mality of the basis functions g: 

I g*(p)d + g(k) = - 2Z)(1 - 2L) - (1 - 2L)(1 - 21) + -(2 -21- 2L)A(2 -21- 2L)\k) 

J x 2 2 
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+ -(p|(l - 2/)(l - 2L) - (1 - 2L)(1 - 21) + -(2 - 2/ - 2L)(A - A f )(2 - 2/ - 2L)|fc) 



-z2vr5(- 
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1) = -i2ir5(r] - f ) . 
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Now for the calculation of the overlap matrices 

e(p,p~,fc,fc - ) = i f d + g*(p)f(k) = i(p\2l-2L + (l-l-L)A*\k) (B.3) 

J x~ ,x 

+ / —{p\d 2 e~ ix ~ [ ~~~\l-2l)\k) 
Jx-<o 2tt 

+ / ^(p| J D 2 (l-2L)e~ i ^ [ ~"~ 1 |fc) 

Jx->0 27T 

= i(p\2l - 2L + (1 - I - L)A^\k) 

+ - p - . (Pl(l " 2Z)|fc> - <<P|(1 " 2L) 1 ■ |fe) . 

By inspection we see that to get the matrix $ we need to substitute k~ — > — /c~ and also A; — > — fc. 
Thus we have 

$(p,p~; — fc, fc _ ) = i(p|2i - 2L + (1 - I - L)A t |/c) 

+ - ' , <P|(1 " 2Z)|fc> - i(p|(l - 2L) 1 . \k) . 

We have been a little cavalier about the non-commuting factors D and d in the above expressions. 
It is clear by construction that all factors of D come from the function g* . Thus all the above 
expressions should be understood as ordered such that all factors D are to the left of factors d. 



C. Appendix - Perturbative expansion. 

Our aim in this appendix is to show explicitly that the 0(Y 2 ) terms discussed in Section 3 have 
their exact counterpart in perturbative calculation. To do this we have to expand our expression for 
the matrix K to first order in j. Our discussion will be for the wave function rather than for the S- 
matrix, but the extension is straightforward. 

We thus consider the expansion of the operator A to O(j). This is straightforward and we give 
some details below. In all the expressions below the field b has to be understood as bf = ^pj a - We 
start by expanding the matrix $ 

HP, m -k, = i(p\2(l — L) + (1 — / — L)At|fc) + y^z^pKi - 2l)\k) - i( P \(l - 2L) — ^_ v \k) . 

W (C.l) 

We have neglected the it terms in this expressions since the poles are never reached for physical 
values of rapidity. 

2(1 - L)f 3 = 2r b (ikjjfl + q^dj - di± {d k , b%] Is,-) . (C.2) 

Further 

4? = f cab dA [d t , bf] ±dj - 2r%±d . (C.3) 



-29- 



And 



(l-l- L)ijA 
Thus the first term in <J> is 



ad a jdb rcab 

ij ^kj ~ J 



2(1 ~ L)f 3 + (1 - I - L)#A{f = r h (id^ + [d k , b%] 



' Q2"3 1 ~ L Q2 J fcJ 02 ^ 

The second term in <I> is of order one, while for the third term we need 

db 



(C.4) 



(C.5) 
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1 + et-v 
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d 2 1 + et-v 
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Thus finally 
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d 2 3 
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(C.6) 
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1 + et-v 
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1 e* - " 



6 ij -2d i -^d j ){d k ,b c k } a2l + e€ _ t? 



(C.7) 



The expansion of matrix G is trivial, since we only need its leading order term 



e 



db 
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ab 



1 - e?-*? - ie 1 - e?"'? + ze 
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Zdi^djS"* ) <M; - ,/) 



(C.8) 



This yields 



A#(x,i7;i/,£) = / er fe ^(x,r^,c)^(z,C,y,£) 



> — lad 



if^W-di^-b^d, 



cab 



(l + e?"'?) 



■^ 5 iMi d h h i}-^\ 



(C.9) 



In this calculation we are only interested in terms which contribute 0(Y 2 ) when A is integrated over 
both rapidities. The second term in eq.( |C.9| ) exponentially decreases for £ — r\ — ► ±oo and thus does 
not contribute. We can then limit ourselves to 



A%(x,r,;y,0 « Ag>(z,y) = -igf cab <j y)6, c (y) + b?(x)-^{x, y)^(y) 



(•cafe 



1 



1 



02 



'02 



(CIO) 
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where we have also restored the coupling constant, omitted from the previous formulae for simplicity. 
We now concentrate on the expansion of the evolved wave function 



1^ = e *J.,.^wfa(*.o^I>,o) e -u^^ . 



(C.ll) 



Expanding in the Fock basis we find 
1 



l-Y- 



|0) 



Y 



i\ bt(x)al a (x,rj)- ig — bt{x)>tf(x,y)a} a (y,Ti) 



|0) 



+ 



X^(x,y)4 u (x, V ) a y(y,0 - \ / 6?(x)65(y)at°(x,T,)at 6 (y,C) 



x,y,r],l 



|0> 



Hence the one gluon component is can be written as 



(a,i,x,r,\V) =ibUx)+g^f abc bHx) J ^{x,y)b%y) = ig^f + g^f^^ Z )j\ Z )^y)f{y) 

(C.12) 

where now the charge j is normalized so that the charge of a single gluon is of order unity. The two 
gluon component of the wave function is 

diTp{x, u)f {u)dj ^ (y, v)j b (v) 



+ 



^9y abc J u ^d~(x,y)d j ^(y,u)j%u) - d j ^(y,x)d l ^(x,u)f(u) >[ j . (C.13) 



These expressions have a very simple diagrammatic interpretation. 

The diagrams contributing to the one gluon state components are depicted in Fig. 4 while to 
the two gluon components in Fig. 5. Clearly the first diagram on Fig. 4 represents the first term in 
eq.( |C.12] ) and the second diagram - the second term in eq.( C.12| ). The second term has an extra 
power of Y due to the additional loop integral. Thus when calculating the 5-matrix, the square of 
the first term is of order 0(Y), while the cross term is of order 0(Y 2 ). This is precisely the term we 
have encountered in our calculation in the body of the paper. 



-f 




Figure 4: Diagrams contributing to the one gluon component of the wave function in perturbation theory. 

Similarly Fig. 5 represents the two terms in eq.( |C.13| ). Both terms here when squared, give 
contribution of 0(Y 2 ) to the S matrix. 
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Figure 5: Diagrams contributing to the two gluon component of the wave function in perturbation theory. 
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